Correlations beyond the horizon 
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We have imaged spontaneously created arrays of vortices (magnetic flux quanta) , generated in a 
superconducting film quenched through its transition temperature at rates around 10 9 K/s. Spon- 
taneous appearance of vortices is predicted by Kibble-Zurek and by Hindmarsh-Rajantie models of 
phase transitions under non-equilibrium conditions. Differentiating between these models requires 
a measurement of the internal correlations within the emerging vortex array. In addition to short 
range correlations predicted by Kibble and Zurek, we found unexpected long range correlations 
which are not described by any of the existing models. 



Any physical system undergoing a phase transition 
within a finite time interval is necessarily driven out of 
equilibrium. One model that describes the dynamics of 
such phase transitions is the Kibble-Zurek scenario. The 
model was first suggested by Kibble [H, Q for cosmolog- 
ical phase transitions which occurred in the early uni- 
verse. Kibble proposed that fast cooldown of the uni- 
verse, through a critical temperature T c , leads to forma- 
tion of initially isolated domains of a new ordered phase. 
The typical size of a domain, £ , inside which the emerg- 
ing order is coherent, is a product of the speed of light 
and time needed to complete the phase transition. Conse- 
quently, different domains are uncorrelated due to causal- 
ity. As these disjoint domains coalesce, the mismatch of 
the order parameter between different regions leads to 
the appearance of topological defects. The domain size 
determines both the typical distance and the correlation 
length between topological defects. Zurek [3] proposed 
terrestrial tests of this model by examining analogous 
situations in condensed matter systems having the same 
symmetry of the order parameter. The analogous "speed 
of light" in condensed matter systems is the velocity of 
propagation of the order-parameter. Since, some aspects 
of the Kibble-Zurek (KZ) model has been tested in a va- 
riety of physical systems, including liquid helium [3, 
liquid crystal s [y| , superconductors [7|, |8[ and Josephson 
junctions [l, fief . Thus, the Kibble-Zurek model is a 
universal theory of defect formation whose applications 
range from phase transitions in grand unified theories of 
high-energy physics to phase transitions observed in dif- 
ferent condensed- matter systems. An alternative mecha- 
nism of spontaneous vortex formation in superconductors 
was proposed by Hindmarsh and Rajantie [jj], [l2[ (HR). 
According to this mechanism, thermal fluctuations of the 
magnetic field freeze inside the superconductor during 
the transition, creating domains of magnetic flux with the 
same polarity and characteristic size. In contrast, the KZ 
model predicts only short range correlations, with neigh- 
boring defects having different topological charge. In a 
superconductor, where the topological defects are vor- 
tices carrying a quantum of magnetic flux, this means 
that adjacent vortices should have a different polarity. 
The presence of spontaneously generated topological de- 



fects was observed in several experiments. Regarding the 
more sensitive testing of correlations predicted by these 
models, there is only one experiment on liquid crystals 
in which an array of topological defects was actually im- 
aged. However, the amount of data was insufficient to 
detect correlations beyond nearest neighbours. The ob- 
jective of our experiment is to decide between these two 
models by imaging spontaneously formed arrays of vor- 
tices in a superconductor and measuring their correla- 
tions. 

Investigating the spatial distribution of the vortex ar- 
ray requires a technique capable of imaging a relatively 
large area with a micrometer resolution. Furthermore, 
the statistical nature of the problem requires accumulat- 
ing many such images. The only practical method to 
achieve both these goals is Magneto-Optical (MO) imag- 
ing. Resolving single vortices however requires push- 
ing the technique to its limits, and until now only one 
group has achieved this goal [13[. We have developed 
a high resolution Magneto-Optical system specifically 
for this experiment. Our system is described in detail 
elsewhere [14]. Briefly, having the magneto-optical indi- 
cator as close as possible to the superconductor's surface 
and a cryogenic design suppressing vibrations allowed us 
to achieve the best spatial resolution (0.8/ira) demon- 
strated so far by this technique [3] • The superconductor 
sample consists of a 200nm thick Niobium film deposited 
on a sapphire substrate. The film was prepared by DC 
magnetron sputtering and it is superconducting below 
8.9K. The film is patterned into small squares of 200/im 
across, to ensure homogeneous illumination by the heat- 
ing laser pulse and avoid thermoelectric currents inside 
the sample. On top of the Nb film, we deposited a 40nm 
layer of EuSe which serves as the Magneto- Optic sensor. 
Due to the huge magneto optical Kerr effect in EuSe the 
polarization plane of linearly polarized light is rotated 
if local magnetic field is present. By mapping this ro- 
tation at each point of the image, we get an image of 
the magnetic field directly above the surface of the su- 
perconductor. The typical lateral size of a vortex in a 
Nb film is about lOOnra [l5[, so the images of individual 
vortices are diffraction limited. In our setup, observation 
of individual vortices is possible over a large field of view 
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(100 x 100/im 2 ) using relatively short integration time 
(10s), allowing us to acquire hundreds of images during 
an experimental run. 

From our previous experiments [7], we know that ex- 
tremely high cooling rates are essential for spontaneous 
generation of a measurable amount of vortices. No less 
important, fast cooling to low temperatures (far below 
T c ) traps the vortices on pinning centers, preventing mu- 
tual annihilation of a vortex and an anti-vortex. High 
cooling rates are achieved in the following way: first, the 
superconducting film is heated above T c by a short laser 
pulse. The 200nm thick film is deposited on a sapphire 
substrate, which is transparent at the wavelength of the 
laser. Hence, only the film heats up, while the 1mm 
thick substrate remains near the base temperature. At 
the end of the heating pulse, the heat from the film es- 
capes via ballistic phonons into the cold substrate, which 
has a thermal mass of ~ 1000 larger than that of the 
film and acts as a heat sink during the cooldown. At 
low temperatures involved (below 14if ), phonon scatter- 
ing is small enough so that the heat transfer from the 
film into the substrate is ballistic. The timescale of the 
heat transfer is much shorter than the length of the laser 
pulse. Therefore, the cooling rate depends only on the 
decay time of the laser pulse. We used pulse shaping 
techniques to change the decay time. In this way two 
cooling rates of 4 • 10 8 and 2 • 10 9 K/s were achieved. 
The cooling rates were measured using GeAu thin film 
resistive bolometer. The conductivity of the GeAu film, 
which is temperature dependent, was measured vs. time 
during the laser pulse. Finally, at these cooling rates, 
the mean inter- vortex spacing predicted by the models is 
larger than our optical resolution, so that we should be 
able to observe the entire vortex array. 

Typical images of spontaneously generated vortex ar- 
rays are shown in Fig.l. During these measurements, 
the system was carefully shielded from external mag- 
netic fields, and the measured asymmetry between the 
density of positive and negative vortices is less than 1%. 
It is clear that the vortex arrays appear random. On 
average, the two cooling rates we could use, 4 • 10 8 and 
2 • 10 9 if/s, produced a density of vortices of 6 • 10 5 cm -2 
and 1.3 • 10 6 cm -2 respectively. Those densities are al- 
most 2 orders of magnitude lower than expected from 
the KZ model. This discrepancy was already noticed in 
our earlier experiments [7]. However, the scaling of the 
density with the cooling rate is consistent with the KZ 
model (proportional to the square root of the cooling 
rate). One potential reason for the low density could be 
the mutual annihilation of nearby vortices having oppo- 
site polarities. The closer vortices are to each other, the 
stronger the attractive force between them. Vortices are 
prevented from mutually annihilating by pinning forces, 
which however do not depend on the distance. Therefore 
a critical distance between vortices with opposite polar- 
ity exists, below which pairs of vortices will overcome the 




FIG. 1. Spontaneously created vortices. Typical images of 
spontaneously created flux lines in a superconductor cooled 
at: a) 2 • 10 9 if/s, and b) 4 • 10 8 K/s. The intensity is pro- 
portional to the local magnetic field. Bright and dark spots 
represent vortices with opposite polarity. 

pinning force, merge and annihilate. We determined this 
critical distance by repeating the experiment under var- 
ious external magnetic fields (Fig. 2). When the density 
of the vortices due to external field is such that the dis- 
tance between them is less than critical, all the vortices 
having an opposite polarity should annihilate. This crit- 
ical length turns to be lpm or less, significantly smaller 
than the scale of short range correlations. Hence, annihi- 
lation should not affect any of the correlations observed 
in our experiment. 

We used images like Fig.l to determine the correla- 
tions between the vortices. In order to increase the sta- 
tistical ensemble, the correlation function was averaged 
over 170 such images, with 12,000 vortices in total. We 
first show the density-density correlation function, irre- 
spective of the vortex polarity. This function is defined 
as D(r — r') =< p(r)p(r') > where p(r) is the local vor- 
tex density. The value of p(r) is taken to be 1 at the 
location of a vortex regardless of its polarity and else- 
where. According to the KZ model, if a vortex is created 
at some vertex between different domains, the probabil- 
ity to find another vortex at the nearest neighbour vertex 
is higher by 33% than the average. Consequently, D(r) 
should show a peak at the characteristic nearest neigh- 
bour distance. The correlation function calculated from 




FIG. 2. Images of magnetic flux above the superconductor 
cooled in the presence of external magnetic field. The exter- 
nal field in (a) is 0.2raT and in (b) is 4 times larger, 0.8raT. 
The cooling rate is 10 9 K/s in both panels. The majority 
of the vortices (dark spots) are associated with the external 
magnetic field. Bright spots are vortices with polarity oppo- 
site to the external magnetic field. Due to annihilation, at 
higher external fields the number of vortices with opposite 
sign is reduced. 

our data is presented in Fig. 3. There is indeed a peak 
at r ^ 3/im, a strong evidence for the short range cor- 
relations predicted by the KZ model. Beyond the peak, 
the correlation function decays towards a constant value. 
The nearest neighbour distance is also the best estima- 
tion for a typical domain size £. 

Next we calculate the correlation function taking 
into account the polarity of each vortex. The vortex- 
vortex correlation function is defined as G(r — r') =< 
n(r)n(r f ) > , with n(r) = 1 at the location of a posi- 
tive vortex, —1 at the location of negative vortex and 
elsewhere. The KZ model predicts that nearest neigh- 
bor vortices should have opposite polarities. This should 
manifest itself as a negative peak in G(r), which ac- 
cording to theory [16] should decay exponentially to zero 
(G(r) oc r 2 exp(— r 2 /£ 2 ). In the KZ model, £ is the only 
length scale, and both the decay length of the correlations 
and the characteristic distance between vortices should 
be = £. In contrast, in the HR model neighboring vor- 
tices should have the same polarity and so the peak in 
the correlation function should be positive. 
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FIG. 3. Density-density correlation function D(r). The error 
bars reflect the statistics. The peak at 3/im is evidence for 
strong short range correlations predicted by the KZ model. 

Fig. 4 shows the vortex-vortex correlation function. 
The correlation function was multiplied by r to empha- 
size long range behavior. The nearest neighbour peak 
is negative, which indicates that the KZ scenario is the 
dominant mechanism of vortex formation. The solid line 
is a fit to the theory of [16[ with £ = 3. It indeed ap- 
pears that the nearest neighbour distance (the peak in 
Fig. 3) and the decay length of the short range correla- 
tion (Fig. 4) are the same. However, at distances beyond 
the peak the theoretical correlation function decays to 
zero, while the experimental correlations do not. Sur- 
prisingly, the data show an oscillatory behavior which 
decays as ex r _a , with a = 1 ± 0.5. This oscillatory 
behavior is well outside the error margins. Such oscil- 
lations, if any, are within the noise in Fig. 3. Hence the 
oscillations in Fig. 4 represent long range correlations in 
the polarity of vortices rather than in their density. We 
found that the wavelength of these oscillations decreases 
weakly with increasing cooling rate as well as with ap- 
plied magnetic field. One possibility is that these long 
range correlations result from some local inhomogeneities 
in the sample. To check this possibility, we repeated the 
experiments at several different locations on the film. We 
found that the correlations were independent of the loca- 
tion. Further, local variations in the local properties such 
as the value of T c could potentially modulate the density 
of vortices which however is constant (fig. 3), but not 
the polarity of vortices which is what we see. 

Long range correlations were predicted in the HR 
model pj], [l2[. However, these correlations should de- 
cay as r -4 , instead of r _1 which we observe and the pre- 
dicted domain size is on a mm scale, 2 orders of magni- 
tude larger than the few fim which we see. In cosmology, 
the size of a domain is determined by the speed of light. 
The analogous speed in superconductors is the propa- 
gation velocity of the order parameter [3]. Close to T c , 
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FIG. 4. First moment of the vortex- vortex correlation func- 
tion G(r). The negative peak at short distance reflects vortex- 
antivortex correlations predicted by KZ model. The blue line 
is a fit to theory of [l6j]. 

this speed is about 10 3 m/sec. This speed determines the 
nearest neighbour distance, the only correlation length in 
the KZ model. However, vortices are also coupled to elec- 
tromagnetic field, which propagates with speed c. Such 
coupling may perhaps be linked with the longer range 
correlations, which are currently a puzzle. 
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